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(54) Improvements in or relating to speech coding 



(57) An improved filtering method for use in an 
enhancement filter in a mixed excitation linear predic- 
tion (MELP) speech coder or a postfilter in a codebook 
excitation linear prediction (CELP) speech coder is dis- 
closed which includes two filters. The first filter (62) has 
a transfer function of 
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where P is the predicted value, a and p are scaling fac- 
tors, z is the inverse off the unit delay operation used in 
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the transform representation of the transfer functions 
and sig-prob is signal probability estimator value and 
the second filter (65) has a transfer function of 1- pz~ 1 * 
sig-prob, where \i = a scaling factor. The sig-prob is the 
signal probability value based on a comparison of 
power of the signals in a current frames to a long term 
estimate off noise power in signal probability estimator 
(63). The sig-prob value is 1 rf the power of the signals 
is greater than the noise power plus 30 dB and the sig- 
prob is zero if the power is less them noise power plus 1 2 
dB. Between these two conditions, sig-prob is (log gain 
- 12 dB - noise gain)/18. 
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Description 

TECHNICAL FIELD OF THE INVENTION 
s This invention relates to speech coding and more particularly to adaptive filtering in tow bit rate speech coding. 
BACKGROUND OF THE INVENTION 

Human speech consists of a stream of acoustic signals with frequencies ranging up to roughly 20 KHz; however, 

10 the band of about 100 Hz to 5 KHz contains the bulk of the acoustic energy. Telephone transmission of human speech 
originally consisted of conversion of the analog acoustic signal stream into an analog voltage signal stream (e.g. ( by 
using a microphone) for transmission and reconversion back to an acoustic signal stream (e.g., by using a loud- 
speaker). The electrical signals would be bandpass filtered to retain only the 300 Hz to 4 KHz frequency band to limit 
bandwidth and avoid low frequency problems. However, the advantages of digital electrical signal transmission has 

is inspired a conversion to digital telephone transmission beginning in the 1960s. Digital telephone signals are typically 
derived from sampling analog signals at 8 KHz and nonlinearly quantizing the samples with 8 bit codes according to the 
li-law (pulse code modulation, or PCM). A clocked digital-to-analog converter and companding amplifier reconstruct an 
analog electrical signal stream from the stream of 8-bit samples. Such signals require transmission rates of 64 Kbps 
(kilobits per second) and this exceeds the former analog signal transmission bandwidth. 

20 The storage of speech information in analog format (fa example, on magnetic tape in a telephone answering 
machine) can likewise be replaced with digital storage. However, the memory demands can become overwhelming: 10 
minutes of 8-bit PCM sampled at 8 KHz would require about 5 MB (megabytes) of storage. 

The demand for lower transmission rates and storage requirements has led to development of compression for 
speech signals. One approach to speech compression models comprises the physiological generation of speech and 

25 thereby reduces the necessary information to be transmitted or stored. In particular, the linear speech production model 
presumes excitation of a variable filter (which roughly represents the vocal tract) by either a pulse train with pitch period 
P (for voiced sounds) or white noise (for unvoiced sounds) followed by amplification to adjust the loudness. The term 
1/A(z) traditionally denotes the z transform of the filter's transfer function. The model produces a stream of sounds sim- 
ply by periodically making a voiced/unvoiced decision adjusting the fitter coefficients and the gain. Generally, see Mar- 

30 kel and Gray, Linear Prediction of Speech (Springer- Verlag 1976). 

To reduce the bit rate, the coefficients for successive frames may be interpolated. However, to improve the sound 
quality, further information may be extracted from the speech, compressed and transmitted or stored. For example, the 
codebook excitation linear prediction (CELP) method first analyzes a speech frame to find A(z) and filter the speech. 
Next, a pitch period determination is made and a comb filter removes this periodicity to yield a noise-looking excitation 

35 signal. Then the excitation signals are encoded in a codebook. Thus CELP transmits the LPC filter coefficients, the 
pitch, and the codebook index of the excitation. 

Another approach is to mix voiced and unvoiced excitations for the LPC filter. For example, McCree, A New LPC 
Vocoder Model for Low Bit Rate Speech Coding, Ph.D. thesis, Georgia Institute of Technology, August 1992, divides the 
excitation frequency range into bands, makes the voiced/unvoiced mixture decision in each band separately, and com- 

40 bines the results for the total excitation. A mixed excitation linear prediction (MELP) coefficient vocoder is described in 
an article by A. McCree, et al. entitled "A Mixed Excitation LPC Vocoder Model for Low Bit Rate Speech Coding", in 
IEEE Trans, on Speech and Audio Proc, Vol. 3, No. 4, July 1995. 

Most low bit rate speech coders employ some form of adaptive spectral enhancement fitter or postf ilter to improve 
the perceived quality of the processed speech signal. For example, in the Mixed Excitation Linear Predictive (MELP) 

45 speech coder in McCree, et al. an adaptive pole/zero enhancement filter based on the LPC spectrum is used. Hie 
adaptive spectral enhancement fffter helps the bandpass filtered speech to match natural speech waveforms in the 
formant region. This adaptive filter described above improves the speech quality for dean input signals, but in the pres- 
ence of acoustic noise this filter may actually degrade performance. The enhancement filter tends to increase the fluc- 
tuations in the power spectrum of the acoustic background noise, causing an unnatural "swirling" effect that can be very 

so annoying to listeners. A similar effect takes place in the postfilter of the CELP speech coder. 

In accordance with one object of the present invention an improvement is provided to this adaptive spectral 
enhancement filter or postfilter in CELP which results in better performance in the presence of acoustic noise while 
maintaining the quality improvement of the existing method for clean speech signals. 

55 SUMMARY OF THE INVENTION 

In accordance with one embodiment of the present invention, a filtering method for improving digitally processed 
speech in low bit rate speech or audio signals is provided wherein the filtering is controlled by linear predictive coeffi- 
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dent parameters and the estimated probability that the input frame is speech rather than background noise. In this way, 
the benefits of tittering are realized for clean speech signals without introducing artifacts to the processed background 
noise. 

5 DESCRIPTION OF THE DRAWINGS 

These and other features of the invention that will be apparent to those skilled in the art from the following detailed 
description of the invention, taken together with the accompanying drawings in which:- 

10 Rg. 1 is a general block diagram of a speech communication system; 

Fig. 2 is a block diagram of the speech analyzer of Rg. 1 ; 

Fig. 3 is a block diagram of a synthesizer; 

75 

Rg. 4a-d illustrates natural speech vs. decaying waveforms where 4a illustrates a first formant of natural speech 
vowel; 4b synthetic exponentially decaying resonance; 4c pole/zero enhancement filter impulse response for this 
resonance; and 4d enhance decaying resonance; 

20 Fig. 5 is a block diagram of the adaptive spectral enhancement according to one embodiment of the present inven- 
tion; and 

Fig. 6 is a flow chart of the signal probability estimator. 

25 DETAILED DESCRIPTION OF PREFERRED EMBODIMENT 

The overall low bit rate speech communication system is illustrated in Rg. 1 where the input speech is sampled by 
an analog to digital (A/D) converter and the parameters are encoded and sent to analyzer 600 ami are sent via the stor- 
age and transmission channel to the synthesizer 500. The decoded signals from the synthesizer 500 are converted 

30 back by the digital to analog converter (DAC) to signals for the speaker. Referring to Rg. 2, there is illustrated some 
blocks of the analyzer. The analog input speech is converted to digital speech at converter 620 and applied to a speech 
analyzer which includes an LPC extractor 602, a pitch period extractor 604. a jitter extractor 606, a voiced/unvoiced mix- 
ture control extractor 608, a gain extractor 610, and an encoder 612 for assembling these five block inputs from 602- 
610 and outputs and clocking them out encoded over a transmission channel. At the synthesizer 500 there is the 

35 decoder 536 which decodes the encoded speech from encoder 612 to provide the LPC parameters, pitch period, mix, 
jitter flags, and gain. 

Referring to Rg. 3 there is illustrated a MELP vocoder according to one embodiment of the present invention and 
described in U.S. Patent application serial no. 08/218,003 ffled March 25. 1994 and similar to that in the above cited 
McCree. et al. article. The synthesizer 500 includes a periodic pulse train generator 502 controlled by a pitch period 

40 input from decoder 536, a pulse train amplifier 504 controlled by a gain input from decoder 536, a pulse jitter generator 
506 controlled by a flag input from jitter output of decoder 536, a pulse filter 508 controlled by five band voiced/unvoiced 
mixture inputs from decoder 536. The synthesizer 500 further includes a white noise generator 512, a gain amplifier 
also controlled by the same gain input, noise filter 518 also controlled by the same five band voiced/unvoiced mixture 
inputs, and an adder 520 to combine the filtered pulse and noise. The adder output is the mixed excitation signal e(n) 

45 which is applied to an adaptive spectral enhancement fitter 530 which adds emphasis to the formants to produce e'(n). 
This output is applied to an LPC synthesis fitter 532 controlled by 1 0 LPC coefficients. The output of this is amplified in 
amplifier 533 with gain from decoder 536 and applied to a pulse dispersion filter 534 to get digital synthetic speech. This 
digitized speech is then converted to analog speech for a loud speaker using a digital to analog converter (DAC) 540. 
In accordance with another embodiment of the present invention, the adder output e(n) is applied to the synthesis filter 

so 532 controlled by 1 0 LPC coefficients and the output of the LPC filter is applied to the adaptive enhancement filter 530 
to add emphasis to the formants to produce e'(n). 

In accordance with one embodiment of the present invention, the present invention enhances the adaptive spectral 
enhancement fitter 530. The adaptive spectral enhancement filter 530 in the MELP coder is a pole/zero fitter based on 
the LPC fOter coefficients. This adaptive filter helps the bandpass filtered synthetic speech to match natural speech 

55 waveforms in the formant regions. Typical formant resonances usually do not completely decay in the time between 
pitch pulses in either natural or synthetic speech, but the synthetic speech waveforms reach a lower valley between the 
peaks than natural speech waveforms da This is probably caused by the inability of the poles in the LPC synthesis fitter 
to reproduce the features of formant resonances in natural human speech. There are two possible reasons for this prob- 
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I em. One cause could be improper LPC pole bandwidth; the synthetic time signal may decay too quickly because the 
LPC pole has a weaker resonance than the true formant Another possible explanation is that the true formant band- 
width may vary somewhat within the pitch period, and the synthetic speech cannot mimic this behavior. 

The adaptive spectral enhancement fitter in the above cited McCree article of July 1 995 provides a simple solution 

s to the problem of matching formant waveforms. An adaptive pole/zero filter is widely used in CELP coders since it is 
intended to reduce quantization noise in between the formant frequencies. See article of Chen, et al. entitled "Real- 
Time Vector APC Speech Coding at 4800 bps with Adaptive Post Filtering", in Proc. IEEE Int. Conf. Accost Speech Sig- 
nal Processing, Dallas 1987, pp. 2185-2188. Also see Campbell, et al. entitled "The DOD 4.8 kps Standard (proposed 
Federal Standard 1016)," in Advances in Speech Coding, Norwell, MA: Kluwer, 1991, pp. 121-133. The poles are gen- 

10 erated by a bandwidth expanded version of the LPC synthesis fitter, with a equal to 0.8. According to the McCree article, 
since this all-pole filter introduces a disturbing lowpass filtering effect by increasing the spectral tilt, a weaker all-zero 
fitter calculated with a equal to 0.5 is used to decrease the tilt of the overall fitter without reducing the formant enhance- 
ment. In addition, a simple first-order FIR filter is used to further reduce the low pass muffling effect. In the mixed exci- 
tation LPC vocoder, reducing quantization noise is not a concern, but the time-domain properties of this filter produce 

75 an effect similar to pitch-synchronous pole bandwidth modulation. As shown in Fig. 4, a simple decaying resonance has 
a less abrupt time-domain attack when this enhancement fitter is applied. This feature allows the LPC vocoder speech 
output to better match the bandpass waveform properties of natural speech in formant regions, and it increases the per- 
ceived quality of the synthetic speech. 

As discussed above, the poles of the enhancement filter are the poles of the LPC filter shifted towards the unit circle 

20 in the z-plane by a factor of 0.8. In accordance with the present invention, since this all-pole fitter by itself introduces a 
muffed characteristic to the processed speech signal, a weaker all-zero filter is used in cascade to compensate for the 
spectral tilt introduced by the poles. In addition, another zero is included in the filter to further reduce spectral tilt Chen, 
et al. in U.S. Patent Na 4,969,192, entitled, "Vector Adaptive Predictive Coder for Speech and Audio," used a second 
fitter in a postfilter in a CELP speech coder. 

25 The problem with this existing method is that it increases fluctuations present in acoustic background noise. The 
new method, taught herein, adapts the strength of the spectral enhancement filter based on an estimate of the proba- 
bility that the current input frame is speech rather than background noise. This probability is estimated by comparing 
the power in the current speech frame to a long-term estimate of the noise power. To prevent possible discontinuities 
from switching the enhancement filter "ON" and "OFP, the strength of the filter gradually varies from no filtering at all 

30 to full spectral enhancement over a range of signal probabilities. 

Referring to Fig. 5 there is illustrated a block diagram of the improved enhancement fitter according to the present 
invention. The mixed excitation signal e(n) is applied to fitter 62 which is controlled by the LPC coefficients and which 
has the transfer function of 



where z is the inverse of unit delay operator z"\ a and p are coefficients empirically determined with some trade-off 
40 between spectral peaks producing chirping and not achieving spectral enhancement. The prediction fitter coefficients 
1-P(z) is equal to the analysis filter coefficients A(z). The frequency response in Hz is the difference between the fre- 
quency responses of two all pole fitter as: 

45 20 log He*)! = 20 log ^--20 log — ~— 

li.PC.-M '"^ 



In the prior McCree article, the values for the enhancement f Bter comprised a first filter, where 0=0.5 and a = 0.8 
so and a second fitter of a transfer function of 1 - pz~ 1 . According to the present invention for the first fitter, the signal prob- 
ability (sig-prob) value from the signal probability estimator 63 is multiplied (*)to the p of 0.5 and multiplied (*) to the a 
of 0.8, or p= 0.5 * sig-prob (signal probability as measured at estimator) and a = 0.8 * sig-prob at the filter 62. The output 
of filter 62 is coupled to a second filter 65 which has the transfer function of 1-jiz" 1 multiplied (*) by sig-prob where ji is 
typically 0.5 multiplied by (*) k(1). The term k(1) is the first reflection coefficient The signal probability estimator 63 is 
55 responsive to the gain from the analyzer (610 in Fig. 4 decoded from 536 of Fig. 2) to determine if the power in the cur- 
rent frame compares to a long term estimate of the noise power. A flow chart of the estimator is shown in Fig. 6. The 
estimator 63 sets some time constants and step sizes and then compares the log of the gain to noise gain +30 dB. If 
the power level is greater than noise gain 430 dB, set sig-prob to 1 and if less than noise gain +1 2 dB, set the sig-prob 
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to zero to have no filtering. In this way, the filter is applied if a signal is present but not if noise is present If the gain is 
between these extremes the sig-prob value is equal to (log-gain - 1 2 dB - noise gain) divided by 1 8. This is a linear ramp 
value of between 0 and 1 between 12 dB and 30 dB. 

This "sig-prob" becomes the multiplier for a, 0 and \i. The time constants are selected to average out the voice sig- 
nal and approximate the value of the noise floor. 

In a real-time implementation of a 2.4 kb/s MELP coder running on a TMS320C31 DSP chip, this improved adaptive 
spectral enhancement method results in a clear improvement in speech quality for noisy input speech, while maintain- 
ing the same quality as the existing method for clean input signals. 

The estimator 63 may be part of the processor chip running code following the pseudo code below: 

* Estimate average noise gain from log gain for current frame 

time constants/step size 
up = 0.0675; 
down = -0.27; 
min = 10; 
max = 80; 

if (log gain > noise gain+up) 

noise gain = noise_gain+up; 

else if (log gain < noise_gain+down) 
noise_gain = noise gain + down; 

else 

noise_gain = log gain; 

/* Constrain total range of noise gain */ 
if (noise gain < min) 

noise gain = min; 
if (noise gain > max) 

noise_gain = max; 

* Estimate current frame signal probability by comparing to noise 
power 

if (log_gain > noise gain + 30dB) 

sig_prob = 1.0; 
else if (log gain < noise gain + 12 dB) 

sig — prob = 0.0; 

else 

sig_prob = (log gain - 12 - noise gain) /18; 
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* Calculate postfilter coefficients 

5 pf_num = bw_expand (lpc_coeff, sig_j?rob*0. 5) ; 

pf_den = bw_expand (lpc_coeff, sig_prob*0. 8) ; 

tilt__cof = [1, -sig_prob*k [1 first reflection 
io coefficient] ] ; 

* Apply adaptive spectral enhancement filter to excitation signal 

15 

filter (excitation, pf_ num, pf_den) ; 
filter (excitation, tilt_cof) ; 

20 We note that this method can easily be applied in other speech coding applications where spectral enhancement 
or postfittering is desired. 

Chen, et al., U. S. Patent Number 4,969, 192 cited above described a post filter where the values for the first filter 
are 0 = 0.5 and a = 0.8 and the second fitter transfer function is 1 -jiz~ 1 . In accordance with the teachings herein the short 
delay post fitter 32a when modified as discussed above to account for the estimated probability is speech rather than 
25 background noise such that for the first fflter p = 0.5* sig-prob and a = 0.8 * sig-prob. The second fitter would have the 
transfer function jiz~ 1 * sig-prob, where \i is 0.5 * k(1) where k(1) is the first reflection coefficient. 

Although the present invention and its advantages have been described in detail, it should be understood that var- 
ious changes, substitutions and alterations can be made herein without departing from the spirit and scope of the inven- 
tion. 

30 

Claims 

A filtering method for improving digitally processed speech signals; 

generating a signal probability estimator value based on a comparison of power of said speech signals in a cur- 
rent frame to a long term estimate of noise power; 

first filtering said speech signals by a delay controlled by linear predictive coefficients and said signal probabil- 
ity value; and 

second filtering said speech signals by a transfer function of the form 1 -jiz" 1 * signal probability value where n 
is a scaling factor and, z~ 1 is a unit delay operator. 

The filtering method of Claim 1 , where said step of generating a signal probability estimator value comprises gen- 
erating a signal probability estimator value is 1 of log gain of said power of said signals is greater than noise power 
plus 30 dB. 

The filtering method of Claim 1 or Claim 2, wherein said step of generating a signal probability estimator value com- 
prises generating a signal probability estimator value if zero if said power is less than noise power plus 12 dB. 

The filtering method of any proceding, wherein said step of generating a signal probability estimator value com- 
prises generating a signal probability estimator value equals to (log gain - 1 2- noise gain)/1 8 if said power is greater 
than noise gain plus 12 dB and less than noise gain plus 30 dB. 

The filtering method of any preceding Claim, wherein said first filtering step comprises a transfer function of 

55 
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where P is the predicted value, a andp are scaling factors, z is the inverse of the unit delay z, and \i is a scaling 
factor. 

6. The filtering method of Claim 5, wherein said first filtering step comprises performing said transfer function where 
a = 0.8, 0 = 0.5. 

7. The filtering method of Claim 5 or Claim 6, wherein said first filtering step comprises performing said transfer func- 
tion where \i is 0.5 * k(1), where k(1) is the first reflection coefficient. 

8. The filtering method for enhancing digitally processed speech or audio signals comprising the steps of: 

buffering said speech or audio signals into frames of vectors, each vector having K successive samples; 
performing analysis of said buffered frames of speech or audio signals in predetermined blocks to compute lin- 
ear predictive coefficients and power in the current frame; 

generating a signal probability estimator value sig-prob based on comparison of the power in the current frame 
to a long term estimate of the noise power; 

first filtering each vector by a delay controlled by said linear predictive coefficient and said signal probability 
estimator value, wherein filtering is accomplished by using a transfer function of the form 



where 1 -P is the LPC coefficient z is the inverse of the unit delay operator used in the transform representation 
of the transfer functions, a and 0 are scaling factors * sig-prob; and 

second filtering by the transfer function of the form 1 - nz~ 1 * sig-prob, where n = scaling factor. 

9. The f fltering method of Claim 8, wherein said step of generating a signal probability estimator value comprises gen- 
erating a signal probability estimator value of 1 if the power is greater than noise gain plus 30 dB. 

10. The filtering method of Claim 8 or 9 wherein said step of generating a signal probability estimator value comprises 
generating a signal probability value of zero if the power is less than noise gain plus 12 dB. 

1 1 . The filtering method of any 8 to 1 0 Claims, wherein said step of generating a signal probability estimator value com- 
prises generating the signal probability estimator value equal to (log gain - 12 - noise gain)/18. 

12. The filtering method of any 8 to 11 Claims, wherein said first filtering step comprises performing said transfer func- 
tion where p is 0.5 and a is 0.8 and in is a5 k(1), where k(1) is the first reflection coefficient 

13. The filtering method of any of 8 to 12 Claims, further comprising setting said sig-prob to 1 if the log gain is greater 
than noise gain plus 30 dB. 

14. The filtering method of any of 8 to 12 Claims, further comprising setting said sig-prob to zero if the log gain is less 
than noise gain plus 12 dB. 

15. A low bit rate speech communication system for transmitting speech signals comprising; 

means for buffering said speech signals into frames of vectors, each vector having successive samples; 
means for performing analysis of said buffered frames of speech or audio signals in predetermined blocks to 
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compute encoded speech including linear predictive coefficients and power in the current frame; 
means for transmitting said encoded speech over a transmission channel, 

a synthesizer coupled to said means for transmitting and responsive to said encoded speech for decoding said 
speech into digital signals; 

a digital to analog converter means responsive to said digital signals from said synthesizer for providing 
speech signals, 

said synthesizer comprising means for enhancing digitally processed speech comprising; 

means for generating a signal probability estimator value sig-prob based on comparison of the power in the 

current frame to a log term estimate of the noise power; 

first filter means for filtering each vector by a delay controlled by said linear predictive coefficient and said sig- 
nal probability estimator value, wherein filtering is accomplished by using a transfer function of the form 



where 1 -P is the LPC coefficient z is the inverse of the unit delay operator used in the transform representation 
of the transfer functions, a and 0 are scaling factors; and 

second filter means for filtering by the transfer function of the form 1 - nz" 1 * sig-prob, where \i = scaling factor. 

16. The system of Claim 15, wherein said signal probability estimator value sig-prob is 1 if the power is greater than 
noise gain plus 30 dB. 

17. The system of Claim 15 or Claim 16, wherein signal probability estimator sig-prob is zero if the power is less than 
noise gain plus 12 dB. 

18. The system of any of Claims 15 to 17, wherein the signal probability value sig-prob is equal to (log gain - 12 - noise 
gain)/18 if the power is greater than noise gain plus 12 dB and less than noise gain plus 30 dB. 

19. The system of any of Claims 16 to 18 wherein 0 is 0.5 and a is 0.8 and n is 0.5 k(1), where k(1) is the first reflection 
coefficient. 

20. The system of any of Claims 15 to 19, wherein said synthesizer includes an LPC filter controlled by LPC coeffi- 
cients. 

21. The system of Claims 20, wherein said means for enhancing is before said LPC filter. 

22. The system of Claim 20. wherein said means for enhancing is after said LPC filter. 

23. The system of any of Claims 15 to 22, wherein said system is a MELP coder. 

24. A fitter for improving digitally processed speech signals; 

means for generating a signal probability estimator value based on a comparison of power of said speech sig- 
nals in a current frame to a long term estimate of noise power; 

a first fitter for filtering said speech signals by a delay controlled by linear predictive coefficients and said signal 
probability estimator value; and 

a second fitter having the transfer function of the form 1 -jiz' 1 * signal probability estimator value where p. is a 
scaling factor and, z" 1 is a unit delay operator. 

25. The fitter of Claim 24, wherein said signal probability value is 1 if log gain of said power of said signals is greater 
than noise power plus 30 dB. 

26. The filter of Claim 24 or Claim 25, wherein said signal probability estimator value is zero if said power is less than 
noise power plus 12 dB. 
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27. The filter of any of Claims 24 to 26, wherein if said power is greater than noise gain plus 1 2 dB and less than noise 
gain plus 30 dB the signal probability value equals (log gain - 12- noise gain)/18. 

28. The filter of any of Claims 24 to 27, wherein said first filter has a transfer function 



10 



where P is the predicted value, a and p are scaling factors, z is the inverse of the unit delay z* 1 , and ^ is a scaling 



factor. 



29. The filter of Claim 28, wherein a = 0.8, p = 0.5. 



15 



30. The filter Claim 28 or Claim 29, wherein n is 0.5 * k(1), where k(1) is the first reflection coefficient 
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